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Abstract. We describe the weak lensing methodology we have applied to multi-colour CFH12k imaging of a homogeneously- 
selected sample of luminous X-ray clusters at z ~ 0.2. The aim of our survey is to understand the variation in cluster structure 
and dark matter profile within rich clusters. The method we describe converts a fully reduced CFH12k image into constraints 
on the cluster mass distribution in two steps: (1) determination of the "true" shape of faint (lensed) galaxies, including object 
detection, point spread function (PSF) determination, galaxy shape measurement with errors; (2) conversion of the faint galaxy 
catalogue into reliable mass constraints using a range of ID and 2D lensing techniques. Mass estimates are derived indepen- 
dently from each of the three images taken in separate filters to quantify the systematic uncertainties. Finally, we compare the 
cluster mass model to the light distribution of cluster members as derived from our imaging data. To illustrate the method, we 
apply it to the well-studied cluster Abell 1689 (z = 0. 184). In this cluster, we detect the gravitational shear signal out to ~ 3 Mpc 
at >3-cr significance. The two-dimensional mass reconstruction has a ~10-<x significance peak centered on the brightest clus- 
ter galaxy. The weak lensing profile is well fitted by a NFW mass profile with Mioo = 14. U" x 10 14 Mq, and c = 3.5*qj 
Of 2 = 0.33), or by a power law profile with q = 0.75 ± 0.07 and 9 E = 14'.'6 ± 073 (x 2 = 0.64). The mass-to-light ratio is found 
to be almost constant with radius with a mean value of M/L R = 150 h (M/L)q. We compare these results to other weak lensing 
analyses of Abell 1689 from the literature and find good agreements in terms of the shear measurement as well as the final mass 
estimate. 

Key words. Gravitational lensing: weak lensing - Galaxies: clusters - Clusters of Galaxies: individual (Abell 1689) 



5h ' 1. Introduction Although the study of a single cluster can be instructive, 

we need to study homogeneous samples of massive clusters in 

Clusters of galaxies are the most massive collapsed structures order to better understand cluster physics, test theoretical pre- 

in the Universe. They are located at the nodes of the filamen- dictions and to constrain the cosmological and physical param- 

tary cosmic web, as mapped by the SDSS and 2dF redshift e ters governing the growth of structure in the Universe. Indeed, 

surveys. These massive systems are the focus of both theoret- c i us t e rs are expected to show some variation in their proper- 

icaWe.g. |Eke et alj[l99g |Bahcall etalj|l997t |Viana & Liddji t i es , in particular in the amount of substructure and their merger 

1998) and observational studies. The aim is to better under- historv. which can be direcflv nrobed bv measuring their mass 



history, which can be directly probed by measuring their mass 
stand cluster formation and evolution and thus it is important to distribution. Thus to obtain a representative view of the proper- 
quantify their physical properties as precisely as possible (e.g. t i es of clusters, a fair and statistically-reliable sample of cluster 
mass distribution, mass density profile, importance of substruc- needs to be studied. 

ture, etc.). Different techniques such as galaxy dynamics, X-ray . , , 

. . , „ „ , , . , „ .... In order to obtain a better understanding of the mass dis- 

emission, the Sunyaev-Zeldovich effect or gravitational lens- ., . , , ° 

, , , , , , . , „ , tnbutions on small and large scales in clusters, we have se- 

mg, are available to probe the physical properties of clusters. ° j ■ 1 

. , , , lected a sample of 11 X-ray luminous clusters (Czoske et al. 

Gravitational lensing is a particularly attractive method as it is i „ i r- — rr^ — nirrrn • . . • . . ~ ', — 

,. t , .. t * * . i. t ., . . . , 2003; Smith et al. 2005) identified i n the XBA Cs sample (X- 

directly sensitive to the total mass distribution irrespective of 1 * 1 *' r i : 1 1 n 

.. , . , . . , , k, „. J lrinri , ray Brightest Abell-type Clusters: Ebeling et al. 1996D. All 

its physical state (see the review by Mellier 1999). J ° Jr 1 — jf* _, 

these clusters have X-ray luminosities of Lx > 8 X 10 erg s 

in the 0. 1-2.4 keV band, and all lie in a narrow redshift slice at 

Send offprint requests to: S. Bardeau z ~ 0.2 (from ZA2218 = 0.171 to ZA1835 = 0.253). As XBACS 
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Fig. 1. The full 42' x 28' CFH12k R-band image of Abell 1689. The thick contours represent the number density of bright 
galaxies selected in the R-band: the first contour corresponds to 5 objects per square arcmin, increased by steps of 5 units. The 
maximal density is 36.5 galaxies arcmin -2 in the cluster center. The thin contours represent the mass density reconstruction with 
LensEnt2 and an ICF of 180" (see Sect. l4.2l for more details). Contour levels are respectively 2, 3, 5, 7 and 9<x, while the peak 
value corresponds to a mass density of 1 100 li7o M0.pc~ 2 . North is to the top, East to the right. 



is restricted to Abell clusters jAbell et all 1 19891) . it is X-ray 
flux-limited but not truly X-ray selected. However, a compar- 
ison with the X-ray selected RO SAT Brightest Cluster Sample 
(BCS: Ebelin getalJl998H2000l) shows that ~ 75% of the BCS 
clusters in the redshift and X-ray luminosity range of our sam- 
ple are in fact Abell clusters. Hence, our XBACs sample is, in 
all practical aspects, indistinguishable from an X-ray selected 
sample. 

Using the CFH12k wide field camera ( Cuillandr e et alJ 

2000) mounted at the Canada-France-Hawaii Telescope 
(CFHT), we imaged all 1 1 clusters in our sample in the B, R 
and I bands. In the present paper we present the weak lensing 
methodology we have applied to analyse these images, using 
our observations of Abell 1689 as a test case. 

The first step of any weak lensing work is to correct the 
observed galaxy ellipticities for any observational smearing: 
circularization or anisotropy due to the point spread function 

(PSF). The classical approach to do this is the so-called KSB 

i — — w- LJ - — I 

method ( Kaiser et al. 1995), implemented in the imcat software 

(see also Luppino & Kaiser 1997i Rhodes et al . 2000; Kaiser 



2000). The basic idea is to relate the "true" ellipticity of the 
background sources to the observed ellipticity through polariz- 
ability tensors, which include the smearing effect of the PSF, 
possibly with anisotropic components. In practice these can be 
computed through the combination of the second order mo- 
ments of the light distribution of the galaxies and the PSF it- 
self. However, in this paper we will use an inverse approach 
through a maximum likelihood or Bayesian estimate of the 
source galaxy shape c onvolved by th e local PSF (this method 
was first proposed by Kuiiken 1999). Both the galaxy shape 
and the local PSF are modeled in terms of sums of ellipti- 
cal Gaussians. This approach is implemen ted in the softwar e 
Im2Shape which has been developed by Brid le et all J200 lit . 
The main advantage of Im2Shape is that it provides estimates of 
the uncertainties of the recovered parameters of the sources and 
these uncertainties can then be included in the mass inversion. 

In the weak lensing limit, the ellipticities of background 
galaxies give an unbiased estimate of the shear field induced 
by the gravitational potential of the foreground cluster. The es- 
timate is inherently noisy due to the shape measurement errors 
and the intrinsic ellipticities of the galaxies. Several methods 
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have been proposed to reconstruct the mass density field (or 
the potential) of the foreground structure from the measured 
shear field. Non-parametric methods are usually best to pro- 
duce a mass-map, necessary to identify mass peaks. They can 
also be used to estimate the cluster mass profil e by means of 
the aperture mas s densitometry method ( Fahl maii et al.lfl994l 
Schneiderl ll996l) . On the other hand parametric methods are 
best to constrain the cluster mass profile and total mass by fit- 
ting a radial shear profile to the galaxy ellipticities. 

To illustrate the various methods and techniques used, 
we apply our procedure to one particulary well-studied clus- 
ter from our sample, Abell 1689. Abell 1689 at z = 0.184 
is one of the richest clusters (R = 4) in the Abell cat- 
alog. Abell 1689 is a powerful cluster lens and has been 



jTvson et alJll990t 


Tvson & Fischer 1995; 


Tavlor et alJ 


1998; 


IClowe & Schneider 


2001:lKimr et al 


||2002 


). It has also been 



Newton JAndersson & Madeiskill2004l) . The central structure 
of this cluster is compl ex: from the redshift distribution of 66 
cluster members Girard fet alJ Jl99 7i find evidence for a super- 
position of several groups along the line of sight to the cluster 



center which explains the extraordinarily high velocity disper- 
jkms -1 . Czoske (2004) has recently obtained 



sion of 23 55 +238 >— 



-183 

a new large dataset of more than 500 cluster galaxy redshifts 
in this cluster, which will help elucidate the galaxy distribution 
Abell 1689. Preliminary analysis of these data shows that the 
large scale distribution of galaxies in and around Abell 1689 
is in fact rather smooth and that significant substructure seems 
confined to the very center of the cluster. Thus, even though 
the cluster has clear substructure, it may still be reasonable to 
model the large-scale mass distribution of the cluster with sim- 
ple models, such as the "universal" mass profile proposed by 
iNavarro et alJdl997l) (NFW). 

This paper is organized as follows: |2]briefly presents the 
observations of Abell 1689 used in this paper and gives a sum- 
mary of the data reduction procedure and the conversion of 
the reduced data into catalogues that can be used in the weak 
lensing analysis. In {j3]we present the measurement of galaxy 
shapes and correction for PSF anisotropy using Im2Shape. 
In §|4] we convert the galaxy shape measurements into two- 
dimensional shear maps and radial shear profiles. explains 
how we model the lensing data using both ID and 2D tech- 
niques. In §|6]we compare the projected mass to the light distri- 
bution in the cluster. Finally in {J7]we discuss our method and 
results. In a separate paper (Bardeau et al. 2005, in prep.) we 
will present a more extensive analysis of the mass distribution 
in Abell 1689 combining weak and strong lensing mass mea- 
surements. 

We assume H = 70kms-' Mpc" 1 , Q ra = 0.3, Q A = 0.7. 
At z = 0.18, 1" corresponds to 3.09 kpc (and 1' to 185 kpc). 



pixel scale of 07205. The log of the observations of Abell 1689 
(aj20oo = 13 h l l ra 30 s , 5j2ooo = -01°20'28") is summarized in the 
first part of Tabled 

2.1. Data reduction 



For a detailed description of the data reduction see ICzoskd 
(2002). Here we just give a brief outline. Pre-reduction of the 
CFH12k data was done in a standard way using the iraf 1 pack- 
age mscred ( Valdes 1998) for bias subtraction and flat-fielding 
using twilight sky images. 

Fringing in the I band images was removed by subtracting 
a correction image constructed from eight science images from 
different fields taken during the same night, after masking any 
objects detected in the images. The appropriate scaling for the 
fringe correction was determined interactively. 

Weak lensing applications demand precise measurements 
of the shapes of faint galaxies and therefore precise relative as- 
trometric alignment of the individual dithered exposures of the 
field (~ 6" in our case). A transformation is needed between 
each chip of the input image and a common astrometric output 
grid which has to account for the position of the chip in the fo- 
cal plane, rotation, variations in the height (and possibly tilt) of 
the chip surface with respect to the focal plane, as well as any 
optical distortion induced by the telescope and camera optics. 
Fourth order polynomials were found to be sufficient to model 
these effects. The method t hat we have de veloped follows the 
approach described bv lKaiseretail ( ll999l) . 

We use Digital Sky Survey (DSS 2 ) images to define the 
external reference frame for observations, but then minimize 
the RMS dispersion of the transformed object coordinates from 
all the exposures rather than the deviations between the trans- 
formed object coordinates from the corresponding DSS coor- 
dinates for each individual exposure. This approach ensures 
optimal relative alignment of the transformed exposures. The 
resulting RMS dispersion of the transformed coordinates is of 
order 0701, corresponding to 0.05 of a CFH12k pixel, for usu- 
ally » 100 objects per chip. 

The input images are resampled onto the output grid with 
pixel size 07205 (the median effective pixel scale of the 
CFH12k camera) using the software swarp (Version 1.21). 
Pixel interpolation uses the lanczos3 kernel which preserves 
object counts, without introducing strong artifacts around im- 
age discontinuities ( Bertin 2001 ). Fields with a large number of 
exposures (> 10) were averaged after rejecting outliers, those 
with fewer exposures median combined. 

The images were photometricall y calibrated on fields of 
standard stars taken fr om the list of lLandoltl J 19921) with ad- 
ditional photometry by Stetson (2000). Atmospheric extinction 
was determined from sequences of science images spanning a 



2. Observations and Cataloging 

We observed Abell 1689 with the CFH12k camera through the 
B, R and I filters (Fig.[^shows the R-band image) between 30 
May and 2 June 2000. The camera consists of 12 CCD chips 
of 2k x 4k pixels with a total field of view of 42' x 28' at a 



1 IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 

2 http : //www-gsss . stsci . edu/dss/dss_home . htm, 
http : / /cadcwww . dao . nr c . ca/dss/ 
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Table 1. Observing log for Abell 1689. We indicate the number of detections in each filter (B, R or I), their number density 
(expressed in arcmin~ 2 , in parenthesis), and the magnitude cuts for galaxy classification. Estimated average redshifts z and 
/? = (Di s /D s ), with their standard deviations, are given for the faint galaxy catalogues (see H3.3l for details). 



Filter 


B 




R 




I 




Date of observation 


May 30/June 2 2000 


May 30/June 2 2000 


May 30/June 2 2000 


Number of exposures 


4 




5 




5 




Exposure time (sec) 


JOUU 




J\J\JV 








Seeing 


0.91" 




0.85" 




0.88" 




Completeness mag 


24.9 




24.3 




22.6 




PSF anisotropy 


0.032 ± 0.012 




0.071 ± 0.019 




0.064 ± 0.028 




Number of Detections 


34669 


(28.6) 


41067 


(33.9) 


28805 


(23.7) 


Stars 


2223 


(1.8) 


3488 


(2.9) 


2397 


(2.0) 


Galaxies 


25823 


(21.3) 


30189 


(24.9) 


21145 


(17.4) 


Others 


6623 


(5.5) 


7390 


(6.1) 


5263 


(4.3) 


Bright galaxies 


B<22.0 1171 


(1.0) 


R<21.1 2166 


(1.8) 


I<19.3 950 


(0.8) 


Faint galaxies 


22.5<B<25.4 20186 


(16.7) 


21.6<R<24.7 22794 


(18.8) 


19.8<I<23.3 14382 


(11.8) 


Other galaxies 


22.0<B<22.5 




21.1<R<21.6 




19.3<I<19.8 






orB>25.4 4466 


(3.7) 


or R>24.7 5229 


(4.3) 


orI>23.3 5813 


(4.8) 


Faint galaxies z 


1.02±0.42 




1.06+0.42 




0.82+0.35 




Faint galaxies ft 


0.70+0.08 




0.69+0.08 




0.65+0.07 





sufficient range in airmass to allow accurate determination of 
the extinction coefficient. 

2.2. Object detection 

With the reduced and calibrated images in hand, the weak 
shear information must be extracted from the photometric cat- 
alogues. The analysis of the images involves a number of steps 
that we describe in detail below. These various steps are con- 
trolled in (as much as possible) an automatic way using dif- 
ferent perl scripts which allow a simple and easy handling of 
catalogues and can easily call external programmes. 

In the present paper we first treat the images taken in the 
three filters B, R and I independently . Differences between 
the results obtained from the three datasets are expected due 
to a number of effects. Different seeing in the images affects 
the accuracy of the measurement of galaxy shapes and hence 
the accuracy of the derived shear fields. Different photometric 
depths of the images will change the number density of faint 
background galaxies and thus again the accuracy of the shear 
measurements. Finally, the images sample different wavebands 
of the observed galaxies, which has an effect on the contrast 
between cluster and background galaxies if these are selected 
based on magnitude alone. This independent approach allows 
us to assess the uncertainties introduced by the listed effects. 
Of course it is desirable to eventually combine the information 
present in the three images in an optimal way so as to arrive at 
definitive measurements of the physical properties of the clus- 
ter. A first attempt at this combination is implemented here but 
will be discussed in more detail in a forthcoming paper. 

The first step is to construct a master photometric cata- 
logue of each individual image. For this purpose and to au- 
tomate the procedure as much as possible we have used sex- 



tractor jBertin & Arnoutsl [l996) in a two-pass mode. A first 
run is made to detect bright objects, with a detection level of 
5cr above the background. The average size (full width at half 
maximum, FWHM) of the point spread function (PSF) is then 
easily determined from the sizes of stars. The saturation level 
of the image is also determined in this run. These parameters 
are then fed into a second sextractor run with a lower detec- 
tion level (1.5 cr with a minimum size of 5 connected pixels 
above the threshold). This second output catalogue corresponds 
to the working catalogue. The total number of objects detected 
in each image is given in Table ^ The photometry was com- 
puted using the mag_auto method of sextractor. 

2.3. Star catalogue 

The second step is to extract a star catalogue from the full cat- 
alogue which will then be used to estimate the local PSF. We 
select stars by a number of criteria. First we locate objects in 
the magnitude - /i max diagram (Fig. [2^) where yu max is the cen- 
tral surface brightness of the objects. Stars, for a given flux, 
have the highest peak surface brightness (provided they do not 
saturate the CCD). Hence they populate the "star"-region of 
Fig-Gk limited to a maximum value of the peak surface bright- 
ness by the saturation of the detector, and to a lower value, 
where galaxies start to overlap the star sequence. 

We use an additional cut in FWHM indicated on Fig . [2b : 
objects with FWHM > seeing + 1 pixel are excluded from the 
star catalogue. Note that very compact objects (in the upper- 
right part of Fig.|2k) correspond to cosmic rays or noise defects 
in the overlapping region between chips. They are rejected and 
are put in the "others" catalogue (see Table^. 

Finally, the star catalogue is cleaned one last time (see 33. 1> 
once the star shapes are adequately measured by im2shape. 
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FWHM (pixels) 

b) 

Fig. 2. a) Magnitude-yU max diagram for all objects detected by 
sextractor in the Abell 1689 R-band image. The points in- 
side the parallelogram correspond to stars, the points below to 
galaxies. Points on the upper left correspond to cosmic rays, 
defects and saturated objects, b) FWHM-magnitude diagram 
for all the objects detected by sextractor in the Abell 1689 
R-band image. The vertical black line indicates the average 
seeing value (4.15 pixels for the Abell 1689 R-band image). 
Stars are excluded from the right hatched part of the diagram 
(> seeing + 1 pixel), and galaxies from the left hatched part 
(< seeing - 0.5). 

2.4. Galaxy catalogues 

The third step in our analysis is to compute the galaxy cata- 
logues that will be used to identify the faint lensed galaxies 
and the bright galaxies that are likely to be part of the cluster 
and which will be used to calculate the cluster luminosity. 

Galaxies are selected from the Magnitude-yu max diagram 
(see Fig.^k). First, as for the stars, saturated galaxies are ex- 
cluded. We checked that none of the brightest galaxies in the 
cluster core are affected by this cut which only affects lower 
redshift galaxies. Furthermore, we applied two additional cuts: 
galaxies must have a sextractor class_star parameter lower 



than 0.8 (this removes faint stars or faint compact galaxies from 
the catalogue), and galaxies cannot be smaller than stars, so we 
exclude all objects with a FWHM smaller than seeing-0.5 pixel. 
This blind cleaning is done in a similar way in all three bands. 
These cuts remove most of the defects in the catalogues. 

The galaxy catalogue is then split into three sub-catalogues, 
defined by their magnitude range: one for the brightest galax- 
ies, dominated by the cluster members, one for the faintest 
galaxies expected to be background sources, and the last 
one for the remaining galaxies (intermediate magnitude range 
galaxies or excluded objects). 

The bright galaxies catalogue is defined with respect to the 
apparent m* of cluster galaxies (see ^6.2l for the estimate of m* 
in each filter). In order to achieve good contrast between cluster 
galaxies and the background field population, while still inte- 
grating a fair fraction of the cluster luminosity function, we de- 
fine the bright galaxy catalogue by selecting galaxies down to 
m* +2 for the B and I-band and m*+3 for the R-band (the deeper 
R-band image allows to have a fainter limit). For Abell 1689, 
these correspond to magnitude limits of B < 22.0, R < 21.1 
and / < 19.3. For illustration, a rough estimate of the field con- 
tamination is given for the Abell 1689 R catalogue: outside a 
radius r = 10' the galaxy density measured in the magnitude 
range R < 21.1 is 1.3 galarcmin 2 while the galaxy density in 
an inner radius r = 5' is 5.5 galarcmin -2 . Therefore with our 
selection criteria the field contamination does not exceed 20 
to 25% of the "bright galaxy" catalogue which will be called 
hereafter the "cluster catalogue". After a uniform correction for 
field contamination, it will be used to measure the cluster lumi- 
nosity and derive a light map, providing simple comparisons of 
the stellar distributions between clusters in our survey. 

Fig.|3]shows the colour (R- 1)- magnitude (/) diagram for 
the galaxies matched in both R and I filters The red sequence 
of cluster ellipticals is well defined. The bright galaxies, as de- 
fined above, are plotted as large symbols. These mainly follow 
the colour-magnitude sequence for early-type galaxies within 
the cluster, which indicates that their identification as members 
is likely to be correct. 

A second catalogue is created for the faint galaxies, with the 
following limits: m* + 2.5 < m < nf + 0.5 for the B and I-band 
catalogue and m* + 3.5 < m < nf +0.5 for the R-band catalogue 
(nf is the completeness magnitude which varies from filter to 
filter, see Table [Q. These catalogues are dominated by faint 
and hence probably distant galaxies and are therefore consid- 
ered as catalogues of background galaxies lensed by the cluster. 
The different cuts were adjusted in order to separate the bright 
(foreground) and faint (background) galaxies as much as pos- 
sible without losing too many galaxies (see Fig.0J. 

3. Galaxy shape measurements 

The shapes of stars detected in the images provide our best es- 
timate of the point spread function (PSF), measuring the re- 
sponse of the entire optical system (atmosphere + telescope op- 
tics) to a point-source. The shape of a star includes an isotropic 
component mainly due to atmospheric seeing, as well as an 
anisotropic component caused, for example, by small irregu- 
larities in the telescope guiding. The isotropic component of 
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Fig. 3. Colour-magnitude diagram for the galaxies detected in 
the R and I filters. Magnitudes are the mag_best measurement 
from sextractor, and colours are computed from magnitudes 
measured in a 3" aperture. Larg (red) points are the R bright 
galaxies (as defined in Sect. 12.41 within 300" of the cluster 
centre. 




20 22 
Magnitude 



Fig. 4. Number counts of galaxies in the Abell 1689 R-band 
image in bins of 0.23 magnitude. The thick line corresponds 
to counts of galaxies across the whole field, the thin line to 
galaxies within 300" from the cluster centre. The greyed area to 
the left (right) shows the magnitude selection for bright (faint) 
galaxies. 



the PSF leads to a circularization of the images of small galax- 
ies and thus reduces the amplitude of the measured shear. The 
anisotropic PSF component introduces a systematic component 
in galaxy ellipticities and thus causes a sp urious shear mea- 
surement if not corrected dKaiser et alJll995l) . The geometric 
distortions of the camera and the corresponding instrumental 
shear are corrected during the data-reduction procedure when 
the image is reconstructed on a linear tangential projection of 
the sky on a plane. 

In the case of Abell 1689, which is representative of the en- 
tire survey, the mean anisotropy of the PSF expressed in terms 



■ / M "',v' "fy V <<; r i ' Ht<m« t:> ; 



1 ' 

- 1 000 



-500 
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1000 
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Fig. 5. Abell 1689 R-band PSF map. The length of the vectors 
is proportional to their ellipticity as indicated by the scale in the 
upper-left corner. The origin of the figure is the cluster centre. 
See details in §3.11 



of ellipticity e = (a - b)/(a + b) is much smaller than 0.15 in 
each filter (see Fig. [5}- 

In order to correct for both the PSF circularization and 
the PSF anisotropy, w e use the im2shape software developed 
by bridle et"afl J200ll) . im2shape implements a Bayesian ap- 
proach to measure the shape of astronomical objects by mod- 
elling them as the sum of elliptical Gaussians, convolved by 
the local PSF, which is also parameterized in terms of ellipti- 
cal Gaussians. The minimization procedure of im2shape esti- 
mates the posterior probability distribution of the image given 
the model and the PSF, and Markov Chain Monte Carlo sam- 
pling gives the most probable value for each parameter, with 
the errors linked to the dispersion of the samples. This ap- 
proach i s a pra ctical implementation of the idea presented by 
Kuiiken ( 1999). im2shape is becoming increasingly popular, 
and has been used in a nu mber of weak lensing app lications 
using different instrume nts jKneib et alJl2003(: ICvpriano et alJ 
l2003llFaure et all2004 . 

A detailed comparison between im2shape and the KSB 
method is discussed by Bridle et al. (in prep.). In the follow- 
ing we describe in detail the procedure we implement to trans- 
form the catalogue data into source ellipticity parameters useful 
for a weak lensing inversion. For simplicity, only one elliptical 
Gaussian is used to describe both the shape of the stars and the 
galaxies. Indeed, as shown in Fig.|6] star profiles are well fitted 
by a single Gaussian. Furthermore, orientation and ellipticity 
(the most useful parameters for the weak lensing analysis) are 
relatively insensitive to the model used to describe luminos- 
ity profiles. The a posteriori justification of the validity of the 
choice is demonstrated by the quality of the weak lensing mea- 
surements. 

3.1. Mapping the PSF distribution over the field 

In a first step, im2shape is used to measure the local PSF by 
estimating the shapes of all the stars in the star catalogue. The 
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Fig. 6. Top: a 16 x 16 pixel image showing the averaged shapes 
of the five nearest stars to an arbitrary position (2000,2000) in 
the R-band image of Abell 1689. The contour levels are from 
0.5 to 4.5 in steps of 1. Bottom: A cut along the x-axis of the 
image above, indicated by small squares, and a Gaussian profile 
fit (obtained by im2shape) shown by the solid line. 



resulting PSF catalogue is then inspected in detail. We first re- 
move objects with ellipticity greater than 0.2 which mainly ap- 
pear to be defects between the chips. A second cleaning pass is 
done to remove stars which are very different from their neigh- 
bours: if they are > 2 <x away from the mean value of the local 
seeing, they are automatically rejected from the PSF catalogue. 
The final cleaned distortion map measured from the stars in the 
field is presented in Fig.|5] 

3.2. Faint galaxy shapes 

In a second step, we linearly interpolate the local PSF at each 
galaxy position by averaging the shapes of the five closest stars 
(Fig.|6j. This number of stars is large enough to locally interpo- 
late the PSF, whereas a much larger number would over-smooth 
the PSF characteristics. The efficiency of the PSF measurement 
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Size (ab, in square pixels) 

Fig. 7. Left: The size distribution (a x b) of the deconvolved 
stars. Their average size is consistent with them being point- 
sources. 




Raw ellipticity Deconvolved ellipticity 



Fig. 8. The ellipticity distribution of the faint galaxies in the 
R-band image of Abell 1689. Left: im2shape catalogue with 
no PSF correction. Right: im2shape catalogue with PSF cor- 
rection. The vertical line indicates the average value of the el- 
lipticity. The effect of circularization on the faint galaxies can 
easily be seen. 

and interpolation can be directly tested on the star catalogues. 
Fig.0shows the resulting distribution of the intrinsic sizes of 
stars after deconvolution with the local PSF. They are intrinsi- 
cally much smaller than 0. 1 th of a pixel. 

im2shape then computes the intrinsic shapes of galaxies by 
convolving a galaxy model with the interpolated local PSF, and 
determine which one is the most likely by minimizing residu- 
als. In the end, im2shape's output gives a most likely model for 
the fitted galaxy characterized by its position, size, ellipticity 
and orientation, and errors on all of these. 

Fig. [8] shows how the galaxy ellipticity distribution alters 
after the Im2shape correction; the effect of PSF circularization 
is evident. 

3.3. Mean redshift of the faint galaxies 

Although the photometric catalogues do not contain redshift 
information on the background sources, we attempt to estimate 
this for the population in a statistical sense. This is necessary 
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as the relative distance of the background population and the 
lensing cluster is of prime importance in the quantitative scal- 
ing of the mass distribution from our weak lensing analysis. 
The critical parameter is the mean value of (3 — D L s/Dqs'- 



1 N 



D hS , 



N ^ D 0SJ 



(1) 



where N is the number of faint galaxies in the catalogue and 
£>ls is the angular diameter distance between the lens and the 
source and Dos between the observer and the source. 

To compute ft we have used a photometric redshift cata- 
logue produced from the Hubble Deep Fields (HDF) North 
and South, observed with the Hubble Space Teles cope (HST) 
jFernandez-Soto et aljfl999t IVanzellaet ai1l200ll) . This cata- 
logue, kindly provided to us by S. Arnouts (priv. comm.), gives 
for each object in the HDF-N/S the apparent magnitudes and 
colo urs as well as mea sured spectroscopic redshift if it ex- 
ists ( Vanzella et al. 2002) or a photometric redshift otherwise. 
Similarly, each galaxy in our three CFH12k catalogues (B, R 
or I) has at least one entry in the corresponding photometric 
catalogues. Depending on the number of available entries for 
each galaxy (1, 2 or 3) an automatic search is done in the full 
HDF catalogues for the ten most similar objects in terms of 
magnitude and colours (correcting for the slight differences be- 
tween the photometric systems of the CFH12k and WFPC2 
cameras). Then the average redshift (photometric or spectro- 
scopic if available) of these 10 objects is assigned to the galaxy. 
When photometric measurements are available in all three fil- 
ters for an object, then this procedure crudely mimics a photo- 
metric redshift estimate, while it is a simple statistical average 
of photometric redshifts at a given magnitude limit otherwise. 
Finally, the mean redshift of each catalogue is computed, as 
well as the mean [}. Their values are listed in Tabled 

4. Shear measurements 

We have now measured the "true" shapes of faint galaxies and 
estimated their mean redshift. The lensing equation for galaxy 
shapes can be written as: 

es + g 



1 + g*es 



(2) 



where ei and es are the complex ellipticities of the image 
and the source; g = y/(l - k) is the reduc ed shear; y is 
the shear vector and k is the c onvergence (e.g. lMeilier| [l999: 
iBartelmann & SchneideJ200ll) . Note that both y and k are pro- 
portional to the distance ratio In the weak regime g <k 1 the 
above equation simplifies to: 



ei = e s + g 



(3) 



Assuming that the faint galaxy population lies at our esti- 
mated mean redshift, and assuming that galaxies have random 
orientations in the source plane, it is easy to see that by locally 
averaging a number of ellipticities we have an unbiased esti- 
mate of the reduced shear, allowing us to directly measure the 
mass distribution: 



The bracket () indicate the average of a quantity near a posi- 
tion. However, because of the random orientation of the galax- 
ies in the source plane, the error in the observed galaxy ellip- 
ticities and thus on the estimated reduced shear will depend on 
the number of galaxies averaged together to measure the shear 
dSchneider et all2000l) : 



(i-lglV, 



(5) 



where N is the number of galaxies used in the averaged. cr fs ~ 
0.33 (see Fig. [8) is the dispersion of the intrinsic ellipticity dis- 
tribution, and the ( 1 - \g\ 2 ) factor is the effect of the shear on this 
dispersion. In the weak lensing regime, g is much smaller than 
1 (our measurements reach 0.1 typically, see Fig. II 1 i . and this 
factor can be neglected. The error on the shear measurement is 
then: 



(6) 



<ei> = <g> 



(4) 



Next we will explore different ways to do this averaging 
and constrain the cluster mass distribution. 

4.1. Building the 2D shear map 

The first and simplest test of the lensing influence of Abell 1689 
is to compute the 2D shear maps. To compute the shear maps 
we average the galaxies in cells using the lensing catalogue 
(PSF-corrected faint galaxy catalogue). The cell size is chosen 
so that each cell contains about 35 galaxies. At the magnitude 
depth of the catalogues (~ 20 galaxies arcmin~ 2 ) this number 
is typically achieved for cells sizes of 80" x 80". Averaging 
such a number of galaxies should mean that the measured aver- 
age ellipticity should be small (below 0.03 from Eq. |5J and its 
orientation random in regions with no shear signal. Near mass 
peaks, we expect to see the shear vectors tangentially aligned 
around the centre of mass. Fig.[9]clearly shows that we detect 
this characteristic lensing signal around the cluster core in the 
R-band catalogue of Abell 1689. The signal traced by the co- 
herent alignment of the "average" galaxy shape is represented 
by vectors whose length is proportional to the ellipticity and 
whose orientation follows the mean orientation of the galaxies 
in each cell. Similar shear maps are seen in the two other bands. 

4.2. Reconstructing the 2D mass map 

We use the lensent2 code ( Marsh all et al.l2002l) to compute the 
2D non-parametric mass map of the cluster. lensent2 imple- 
ments an entropy-regularized maximum-likelihood technique 
to derive the mass distribution within the field on a grid. 
The technique consists of a Bayesian deconvolution process: 
a trial mass distribution E(0) is used to generate a predicted re- 
duced shear field through th e convolution of the surf ace mass 
density by a kernel (KS93: iKaiser & Sauireslll993i) . In con- 
trast to KS93, lensent2 cannot produce negative feature in the 
mass maps leading to more physical solutions than can be ob- 
tained from direct reconstructions of the gravitational poten- 
tial if/. Moreover, lensent2 can include information not only 
from the mean shear field, but from each individual lensed 
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Fig. 9. The shear map derived from the R-band image of 
Abell 1689. The cluster centre is marked by a The cir- 
cles have radii of 200"(~ 620 kpc), 780"(~ 2.4 Mpc) and 
1200"(~ 3.7 Mpc) respectively. The inner circle corresponds to 
the region of strong lensing, the second one to the largest circle 
that lies entirely within the CFH12k field, and the outer circle 
marks the limit where the area outside the field becomes signif- 
icant. The shear vectors are computed in cells of 80" x 80", and 
have been smoothed by a Gaussian of 30" width (see details in 



galaxy with its redshift (if known). As clusters of galaxies 
have smooth and extended mass distributions, the values of 
E on the field are expected to be correlated through a kernel 
called the Intrinsic Correlation Function (ICF). In our analy- 
sis, we provide lensent2 with a position, elliptical shape pa- 
rameters (with errors) and an estimate of the redshift for each 
lensed galaxy (we use the mean redshift as explained in ^3.31 . 
There is then only one free parameter in the proceedure, the 
Intrinsic Correlation Function (ICF) which measures the corre- 
lation between mass clumps. We choose a Gaussian ICF, and 
let its width vary. The ICF size is optimized so that the recon- 
structed mass map does not contain a large number of insignif- 
icant small-scale fluctuations, although small ICFs best fit the 
mass peak of the cluster, while large ones best fit the wings of 
the extended profiles. This optimization is performed by maxi- 
mizing the evidence value of each reconstruction, which is the 
probability to observe these data for a given ICF width. For 
more details on lensent2 see lMarshalletal"lll2002h . 

The main cluster mass clump is very well detected by 
lensent2. The code estimates the central surface mass density 
of the peak, and gives its spatial configuration. Note that large 
ICFs smooth the main peak. Reconstructions are computed for 
a large set of ICFs (with scales from 60" to 240"), and the best 
ICF width is found to be near 160-180" for our dataset. An 
illustrative example is shown in Fig.^Swhere the peak of the 
surface mass density is at a value of 1250 h 7 o MQ.pc~ 2 in the 
adopted cosmology, although typical values of the critical sur- 
face mass density for massive clusters at zl ~ 0.2 are roughly 
around E c = 3200 li7o M0.pc -2 for sources at zs ~ 1-0. This is 
because the ICF width used here (180") is much larger than the 
Einstein radius of the cluster (~ 40"). Therefore the smooth- 
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Fig. 10. LensEnt2 mass reconstructions for Abell 1689 from 
the B (top), R (middle) and I (bottom) catalogues. The ICF is 
Gaussian with a width of 180". The cluster peaks are at 1320, 
1250 and 1090 h 70 M .pc~ 2 for B, R and I respectively. White 
(black) in the gray scale is set to lcr (5cr), and contours are at 
2, 3, 4 and 5 cr. cr values are estimated as explained in 34.21 
The scale is in arcsec relative to the cluster centre. A possible 
secondary peak is visible at [-200", -150"]. 



ing process strongly attenuates the central peak density which 
in the case of Abell 1689 is clearly over-critical. 

To assess the significance of the other mass density peaks 
detected in each image we randomize the orientation of the 
faint galaxies in the lensing catalogue, while keeping their po- 
sitions and axial ratios fixed. We perform mass reconstructions 
of 200 randomized catalogues, and in each identify the 15 high- 
est significance mass peaks. The statistics of these 3000 values 
gives a mean noise peak of 1 16 h7o M0.pc~ 2 (99, 85) above the 
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background level (set at 100 1170 Mq .pc 2 in input of lensent2) 
respectively in the R (B, I) images. This value is considered 
as the average fluctuation of the noise peaks, cr. With this def- 
inition, the cluster mass peak is detected at nearly 10cr above 
the background. To be formally correct, prior to randomizing 
their orientations we should also "unlens" the galaxies using 
the shear determined above and applying Eq|2] This has not 
been done yet for simplicity and will be explored in more de- 
tail in the next paper (Bardeau et al. in prep, paper II). However, 
as the lensing induced distortion is almost always very small 
compared to the width of the ellipticity distribution, we do not 
expect that this simplification will affect the estimated signifi- 
cance of the mass peak. 

lensent2 mass reconstructions give many low significance 
mass peaks. For example, Fig.[H)]shows that four clumps reach 
the 2-cr level, but only one is above the 3-cr level (excluding 
the main cluster peak). To check their reality, we can com- 
pare the reconstructions derived from the three filters (B, R 
and I, Fig. I10> . The regions where a mass clump is detected 
in all three images are considered as "real" ones and can be 
compared to the number density map of bright cluster galax- 
ies. Another test is to compare these clumps with any en- 
hancement of the light distribution ( 36. 1> . provided that the 
mass clumps are not associated with "dark clumps". The multi- 
colour approach in our weak lensing survey provides a power- 
ful tool to eliminate most of the inconsistencies created in the 
mass reconstructions from defects in the lensing catalogues. 
For Abell 1689, apart from the mass peak associated with the 
cluster, no other > 3<x peaks were detected in all three filters. A 
possible 2-3<r peak is located 5'South-East of the cluster but no 
obvious counterpart in the galaxy distribution can be identified. 



4.3. The radial shear profile 

We have demonstrated that only one significant mass peak is 
detected in the field of Abell 1689, and that it aligns well 
with the cluster centre indicating it corresponds to the poten- 
tial well of the cluster. In order to quantify the mass of this 
clump we analyse the radial distribution of the shear around 
this peak. Tangential and radial shear components are com- 
puted as a function of the distance to the cluster centre. They 
are averaged in annuli of width AR = R2 - R\ for a mean ra- 
dius /? = (/?!+ R 2 ) 12. AR is kept constant so the S/N of the 
shear roughly decreases as 1 / V/?, in order to keep enough inde- 
pendent points at large radii (a constant S/N requires too large 
annuli at these radii). A quasi-continuous profile is built by us- 
ing a "sliding window" with steps Ar much smaller than AR. 
In practice, we chose AR = 160" (and Ar = 10") for the 
Abell 1689 R image, so about 10 independent points are in- 
cluded in the profile. 

Fig.[^shows the tangential and radial shear profiles from 
the three images of the field. The radial shear should be zero in 
the case of perfect data and a well chosen centre for the annuli. 
In practice, it can be considered as an independent estimator 
of measurement errors (this is also referred to as the 45 degree 
test). In the case of Abell 1689, the radial shear is always lower 



than the tangential shear out to ~ 1 100", arguing for good data 
quality in all three bands. 

Note that in the very centre (R < 70") the shear profile 
appears to drop. The error bars are large due to the low num- 
ber statistics: the area considered is small, reduced still fur- 
ther by the masking effect of the bright galaxies. Moreover 
the depletion of the number density of background gal axies 
in the center due to the magnification bias (Tav lor et al.l l 1998) 
also decreases the number of observable galaxies, although this 
effect is only important in the inner-most annuli. These low 
number statistics does not completely explain the weakness of 
the shear: it can also be under-estimated if unlensed galax- 
ies (such as cluster members) are included in the catalogues, 
which should be more likely towards the cluster core. As a 
consequence, the points inside R = 70" will not be used in 
the modeling of the shear profile. The measurements done by 
Clowe & Schneider ( 2001) using R band images from the ESO 
Wide Field Imager (WFI) are also presented in Fig. [^f or com- 
parison. Our meas urements are quantitatively in good agree- 
ment with those of Clowe & Schneider (2001). Moreover, our 
error bars are smaller and our points less scattered, even if we 
consider the different binnings. This strongly suggests that the 
use of im2shape in the analysis process improves significantly 
the shear measurements. This will be quantified in a forthcom- 
ing paper (Bridle et al. 2005, in preparation). 



5. Modeling the lensing data 

5.1. Description of the mass models 

Three families of mass models are used to fit the measured 
shear profile: a singular isothermal sphere profile (SIS), a 
power l aw profile (Pow) an d finally the "universal" NFW 
profile ( Navar ro et al 1 ll997h . In addition we implemented 
the Aperture Mass Densitometry method (AMD) to compute 
a non-parametric m ass profile from the shear profile itself 
(Fahlman et al. 1994). We recall briefly the basic equations for 
the mass density (p), shear (y) and convergence (k) profiles for 
the three models. 



5.1 .1 . The Singular Isothermal Sphere model 

This is the simplest mass profile used in lensing inversion. It is 
essentially given by the following equations: 



P(r) = 



2nGr 2 



<6) = y{9) = g 

_ Ano 2 D LS 
c 1 D s 



(7) 
(8) 
(9) 



where cr is the velocity dispersion of the cluster. Note that once 
the cluster centre is fixed, this profile depends has only one 
free parameter (6e or equivalently cr), so only one degree of 
freedom is available in the fits. 
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Fig. 11. The tangential shear profile for Abell 1689, in the B, R and I bands. The bin width (AR) is 160". A series of uncorrelated 
points with error bars is displayed for the R band (solid squares). Absolute values of the associated radial component | < e ± > \ 
are indicated in grey at the bottom of the plot, sho wing that the signal is well detected out to ~ 1000" (~ 3 Mpc) from the centre. 
The measurements of lClowe & Sc hneider (2001) are also shown for comparison (open squares). 



5.1 .2. The Power Law model 



be expressed as 



The Power Law model is a generalization of the SIS model, 
where the slope of the mass density profile is a free parameter 
( Schneide ret alJ2000l) . 



7(0) = 

K(0) = 



1 

2-q 



(10) 



(11) 



where q is the slope of the Power Law (q — 1 for the SIS 
model). Once the cluster centre is fixed, this model provides 
two degrees of freedom for fitting. 

5.1.3. The NFW profile 

The NFW profile is derived from fitting the density pro file of 
numerical si mulations of cold dark matter halos (Na varro et alJ 
1 19951 Il997l) . This theoretically-motivated profile is becom- 
ing increasingly po pular in weak lensing analyses of clusters 
jKneib et alJl2003|) as it appears to give a reasonable descrip- 
tion of the observed shear profiles. The mass density profile can 



P(r) 



where 5 C 



and p c = 



(r/r s )(l+r/r s ) 2 



200 



3 ln(l + c) - c/(l 

3# 2 (z) 
8ttG ' 



c) 



(12) 
(13) 
(14) 



r s is the scale radius, H(z) the Hubble parameter and c - r2oo/r s 
the concentration parameter which relates the scale radius to 
the virial radius ^oo- This density profile is shallower than the 
SIS near the center but steeper in the outer parts. Similarly as 
the power law model, once the centre is fixed, it has two degrees 
of freedom: M200 for the normalization of the mass and r s for 
the scale radius, or equivalently 7*200 an d c - The details of the 
analytic expressions for the shear and c onvergence of the NFW 
profile can be found in King et al. (20Q2|). 
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Table 2. Best fit results for the Abell 1689 R-band shear profile. For the SIS, the results are given in terms of Einstein radius 
(9e) and velocity dispersion <t\ os . For the Power Law, 8e is again the Einstein radius and q the logarithmic slope. Finally for the 
universal NFW profile, c is the concentration parameter and r2oo the virial radius. Mjqo is the 2D-projected mass ins ide r^ao in 
units of 10 1 2 h~\ Mq and 6 E is the derived Einstein Radius, (a) refers to the fit results from lClowe & Schneider! ( BOO lb . (b) from 
iKing etal.N2002h . The numbers in italics assume z s = 1.06. 



SIS 


cr lD (km s- [ ) 


e E (") 


X 2 




998 ± 68 


22.4 ± 3.0 


1.98(1) 


(a) 


1028 ± 35 


23.8 ±1.6 




Pow 


q 


0eO 


x 2 




0.75 ± 0.07 


14.6 ±0.3 


0.637 (2) 


(b) 


0.88 


18.0 





NFW 


c 


>200 (/!7<f ' Mpc) 


M 200 (10 12 Mq) 


9e(") 


2 

X~ 




3 5 +05 

J ' -0.3 


1.99 ±0.25 


1410+ 630 
itiu 47Q 




0.334 (2) 


(a) 


6.0 


1.83 


1030 


9.7 




(b) 


4.8 


1.84 


1070 


5.3 





5.2. Weak lensing fit 

Each of the three models presented above is fitted to the data 
with a least square minimization over the parameter space of 
the models. The^ 2 value is then: 



X 



N 



1 N I 

k=l x 



' 4 _ gmodelfe) 



CTk 



(15) 



where N is the number of data bins, and cr^ is the error on the 
tangential ellipticity. The error is computed in each bin as the 
mean error on the tangential ellipticity (e'), weighted by the 
number iVj of galaxies in the bin used to do the measurement: 
cr k = <ov> t / VA^. 

The data in the outer regions at r > r max , where the annuli 
reach the borders of the field, are excluded. In practice, only 
the area where the tangential shear is greater than radial shear 
is included in the fits. Furthermore, as explained in 34.31 we 
also exclude the central part of the cluster. In the case of the R- 
band observations of Abell 1689, the fitted range corresponds 
to radii from r m ; n = 70" to r max = 1 100". 

Table |2] summarizes the results of the fits, and Fig.ll2h dis- 
plays the resulting best-fit models. The lower quality of the fit 
by the SIS profile is easy to understand as it depends on one pa- 
rameter only, contrary to the other two which are represented 
by two parameters. Moreover, the value of the Einstein radius 
deduced from the fit is significantly lower than that measured 
from strong lensing (wh i ch is estimated to be 0e — 41"). Note 
that lClowe & S chneider J200ll) deduced from their weak lens- 
ing analysis a value for the Einstein radius similar to our esti- 
mate. 

The fit with a power law is slightly better than the SIS as 
the slope of the profile appears to be shallower than isother- 
mal, but the Einstein radius is again some 25% lower than ex- 
pected dKin g et alJ |2002) found similar results with an even 
lower Einstein radius). 

The universal NFW profile provides the best fit to our shear 
profile. The concen tration parameter (c) i s slightly smaller than 
the values found bv lClowe & Sch neider (2001) and lKing et al. 



J2002I) . whereas the virial radius r2oo is very similar. The de- 
rived Einstein radius is however quite small and thus this model 
is not a good fit of the central parts of the cluster. 

We conclude that total mass profile of this cluster across 
all angular scales is not well described by any of these simple 
fitting formulae and requires a more complex model, perhaps 
including contributions from the cluster galaxy halos and pos- 
sibly a steeper central mass distribution. 

Fig. 112b shows the projected mass profiles from the previ- 
ous fits computed with the following equation 



M(r) = nr 2 S c Ur) 



(16) 



where k(/) is the mean dimensionless surface mass density in- 
side radius r. 



5.3. The Aperture Mass Densitometry method 

Instead of fitting analytical formulae, we can directly integrate 
the measured reduced shear to determine the relative mass 
profile w ithin the cluster. This direct method has been devel- 
oped by Fahl man et al 1 (119941) and is called 'Aperture Mass 
Densitometry" (AMD). The function f(ri, r 2 ) is defined as the 
difference between the average convergences (or mean pro- 
jected mass densities) inside the radius r\ and within the an- 
nulus between r x and r 2 : 



((ri, r 2 ) = K(r<n)- K(r { < r < r 2 ) 
2 



r yt 

Jn 



lilnr 



The reconstructed mass inside the radius r < r max is then 
M ( (r) = nr 2 S c f (r, r max ) , 



(17) 
(18) 

(19) 



where r max is the maximum radius for which we can measure 
the shear or the radial limit of the data. In the case of our ob- 
servations of Abell 1689, we choose r max ~ 1100", the maxi- 
mum radius where annuli lie entirely within the field of view. 
Regarding Eq.^J M^(r) is only a lower limit to the true mass 
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Fig. 12. Top: Best fitting parameters for SIS, Power law and NFW models, for the Abell 1689 R-band shear profile. One series 
of uncorrelated points is shown (bin width = 160")- Bottom: Deduced mass profiles from these models. The mass profile from 
the Aperture Mass Densitometry Method is also displayed, with a reference radius of 1 100". See text for details. 
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M(r) = nr 2 2 C k{r) and should not be considered as an absolute 
mass determination. 

The AMD mass profile is shown in Fig.^]with the mass 
profiles derived by fitting the various analytical models. As ex- 
pected, we find that the mass estimated from AMD is always 
lower than the parametric mass estimates. 

6. Light distribution and mass-to-light ratio 

6.1. 2D light distribution 

The catalogue of "bright" galaxies is expected to be dominated 
by cluster members, although it may also contain other bright 
galaxies within the field of view. Thus a density map (light den- 
sity or number density) derived from this catalogue can trace 
the morphology of the cluster and any associated structures in 
its vicinity. In the case of Abell 1689, no galaxy over-densities, 
other than the main cluster component, are associated with any 
prominent peaks in the lensing mass distribution (Fig.Q. 

We therefore focus on the distribution of light around the 
cluster centre assuming that the observed over-density is due 
to cluster members. First in order to build a quantitative light 
density map or its radial profile, it is necessary to statistically 
correct the catalogue for the field contamination. Fortunately, 
the CFH12k images are large enough so that at radii beyond 
600" from the cluster centre (2 Mpc at the cluster redshift) 
we can assume that the galaxy density is close to the "field" 
density. The mean number and light densities are therefore cor- 
rected by subtracting their minimal values estimated in the area 
600" < R < 1200". 

Furthermore in order to estimate the total luminosity of the 
cluster and its radial profile, it is necessary to correct for the 
magnitude limit of the catalogue, corresponding to a cut in the 
cluster luminosity function (LF). The incompleteness factor C 
is estimated as follows, the cluster LF is assumed to follo w the 
standard Schechter luminosity function ( Schecht erll976l) : 



dN 

dL L 



-L/LT 



(20) 



Therefore the luminosity integrated in the catalogue down to a 
luminosity L; n f is 



r*+QQ 

= I L 0(L) dL 

•^Anf 

= 0*L* [r(2-a)-r(2-a,L inf /L*)] 



(21) 
(22) 



so the fraction of the luminosity not taken into account when in- 
tegrating within the magnitude limits of the catalogues is writ- 
ten as 



C 



r(2- a, L inf /L*) 
T(2-a) 



(23) 



and the total luminosity is L tot = L cat /(1 - C). 

For the three bands used in this study, we need to estimate 
the two main parameters of the Schechter luminosity function a 
and U . These parameters depend on the choice of filters, on the 
mix of galaxy types, and on the cosmological model. The best 
multi-colour luminosity function determinations are presently 
those from the Sloan Digital Sky Survey (SDSS) early release 



Table 3. The photometric parameters of the luminosity func- 
tion in B, R and I filters for our adopted cosmology: Ho = 
70km s _1 Mpc 1 , Q m = 0.3 and Qa = 0.7. The distance mod- 
ulus is m — M — 39.70 or equivalently the luminosity distance 
is L\ = 872 Mpc, at ZA1689 = 0.18. 1/(1 - C) is the correction 
factor applied to the integrated luminosity of the catalogues to 
determine the total luminosity of the cluster. 





B 


R 


I 


k-correction 


1.06 


0.16 


0.16 


M" - 5 log /770 


-20.47 


-21.83 


-21.54 


a 


1.30 


1.20 


1.25 


m* (z = 0.18) 


20.29 


18.03 


17.32 


1/(1 -C) 


1.28 


1.11 


1.27 



data (Blan ton et alJl200lh . although they correspond to a field 
LF. The SDSS photometric system (u,g,r,i,z) is transformed 
to the CFH12 k (Johnson) sys tem by applying the transforma- 
tions of Fuku gitaet al.Nl996l) . I n this paper we use th e param- 
eters of the LF summarized in Ide Lapparentl J2003I) and ap- 
plied to a Sbc galaxy. Therefore the absolute magnitude M* 
in the R filter is -21.83 in the adopted cosmology and the 
slope is qtr = 1.20. This includes also the k-correction at 
redshift 0.18, computed with the galaxy evolutionary code by 
iBruzual & Charlotill2003l) . 

Finally, the correction factors 1 /( 1 - C) are applied to L cat 
to obtain the total integrated magnitude for the B, R and I cat- 
alogues, with the magnitude ranges defined in Sect. 12.41 The 
precise values derived are summarized in Table|3] 

6.2. Comparison of mass and light: M/L radial profile 

Using the best-fitting NFW model for the observed shear pro- 
file, the M/L profile is computed by dividing the luminos- 
ity profile, estimated from the bright galaxies catalogue, by 
the mass profile. The field contamination in this catalogue is 
estimated by measuring the minimum of the surface bright- 
ness density between 600" and 1200" from the cluster center. 
Fig. I13h displays this integrated luminosity profile for the R- 
band image. Note that the Lr values adopt the correction factor 
discussed in the previous section. 

Fig. 113b displays the Mr(< t)/Lr(< r) profile with error 
bars estimated from the errors on the mass profile only. The 
M/L increases from a low value (near 100 + 101ijo(M/L)q at 
400 kpc from the centre) and flattens out beyond ~ 1 Mpc at 
a value near 160 ± 40 hjo (M/L)q. This behaviour is indepen- 
dent of the filter considered. It does however depend slightly on 
the background correction at large radius, and on the detailed 
mass modeling in the inner part of the cluster. In particular, as 
we derived a relatively small Einstein radius compared to that 
determined from strong lensing, we are likely to be underesti- 
mating the mass in the central regions, which would suggest an 
even flatter M/L profile towards the centre. 

Beyond ~ 1 Mpc the M/L ratio found in Abell 1689 is con- 
sistent with being cons tant with rad ius. This result is similar 
to the findings of Knei b et al.l (|2P03) in their lensing analysis 
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Fig. 13. Left: luminosity profile for the bright galaxy catalogue, for R-band image of Abell 1689, corrected for background 
contamination. See Sect. 16. II for details. Right: Mr/Lr ratio as a function of radial distance from the cluster centre. The mass 
profile is estimated from the best fit NFW parameters. The filled region indicates the errors on the profile. 



of the cluster CI 0024+1654 (z = 0.39), both in the radial dis- 
tribution and in the normalization. For comparison, (M/L)r at 
large radii in the Coma cl uster is found to be 170 + 50 (M/L)p 
from dynamical analysis JGeller et alJl999HRines et alJ20 01). 
Similar profiles for mass and light on 1-5 Mpc scales are 
expected if cluster assembly is largely governed by infalling 
groups and if no strong mass segregation occurs in the cluster. 

In their sample of 12 distant clusters (0.17 < z < 
0.56) ISmail et alJ i 19971) found a mean value of (M/L)f = 
I26+77 7 (M/L)q (h = 0.7) in the cluster cores, where the su- 
perscript all refers to the entire population of the clusters, not 
only early-type galaxies. Given the colour index (V - R) of a 
mean Sa galaxy at redshift 0.18, this corresponds to (M/L) R n = 
102^^ (M/L)q. Since our bright galaxies catalogue is dom- 
inated by elliptical galaxies (FigO, we expect to find a lower 
luminosity thus their M/L value is consistent with our findings. 

7. Discussion and Conclusion 

In this paper, we describe the methodology used to analyze a 
multi-colour wide-field imaging survey of 1 1 X-ray luminous 
clusters. The goal of our survey is to constrain the mass distri- 
bution in clusters of galaxies using weak gravitational lensing. 
The main elements of the data analysis are: the use SExtractor 
for object detection and photometry to provide well-defined ob- 
ject catalogues. A "stars" catalogue is used to determine the 
PSF locally, a "bright galaxies" catalogue is defined to trace the 
distribution of cluster members and a "faint galaxies" catalogue 
is constructed which should comprise background galaxies. 
The magnitude limits of each catalogue are determined with re- 
spect to the observational constraints such as the limiting mag- 
nitudes of the available images as well as physical constraints 
related to the magnitude distribution in the clusters at a given 



redshift. In order to determine the "true" PSF-deconvolved 
shape properties of the background (lensed) galaxies we use 
the im2shape package developed recently for the purpose of im- 
proving the quality of shear measurem ents, including a co rrect 
treatment of the measurement errors jBridle et alJl200lh . We 
then reconstruct the mass distribution by computing the shear 
profile and either fitting it with parametric mass models like the 
NFW mass profile or deducing the relative profile directly with 
the non-parametric Aperture Mass Densitometry method. Both 
methods are found to be consistent. We also sho w a 2D mass re- 
construction using the lensent2 software ( Mars hall et alJ2002l) 
and applying it to the three images taken through the three fil- 
ters. Finally we compute the M/L ratio as a function of radius, 
again in the three photometric bands. The three filters are used 
independently for most of the processing steps in order to con- 
firm the significance of the results (comparison of shear pro- 
files and mass maps). They give quantitatively consistent re- 
sults, further demonstrating the robustness of our method. The 
images in the three filters are used jointly to estimate the back- 
ground galaxies' redshift distribution and so provide a correct 
normalization of the mass determination. 

We apply this method to the well-known cluster Abell 1689 
as a test-case. We find only one significant mass peak in the 
mass reconstructions, corresponding to the cluster itself. This 
is consistent with preliminary results from a large spectro - 
scopic survey of Abell 1689 and its outskirts JCzoskel l2004). 
which shows that the environment of this cluster is remark- 
ably smooth a nd quiet. We also compare our results to pre- 
vious work by Idowe & Schneider! d20Qll) who used an i nde- 
pend ent data set a nd the methods from Kai ser et alJ ([1995) and 
iKaiser & Sauiresl i 19931) for their galaxy shape measurements 
and mass reconstruction. Within the errors both reconstructions 
agree very well. The same is true for the M/L determination, 
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which is consistent with previous findings. Moreover we are 
able to build a M/L profile which in the case of Abell 1689 
shows a near constant behaviour at large radius with a possible 
decrease close to the center. This suggests that mass traces light 
at least in the outskirts of the cluster. The drop of Mj L in the 
cluster centre may be due to an underestimate of the mass in 
the centre, due to increasing contamination of the background 
galaxy catalogue by cluster members diluting the lensing sig- 
nal. The flat M/L profile in the infall region of the cluster in- 
dicates that the association between mass and light has already 
been achieved outside the cluster and the effect of the cluster 
environment on the mass-to-light ratio of infalling galaxies and 
groups is minor. This supports the picture of a hierarchical as- 
sembly of clusters. 

For the results presented here we did not make use of the 
colour information available from multi-band imaging to sepa- 
rate cluster from background gala xies which makes ou r result s 
directly comparable to those of Iciowe & Schneider! d200lh . 
However Clowe (2003) presented an updated mass reconstruc- 
tion for Abell 1689, this time using colours derived from our 
CFH12k images. The colour information resulted in an im- 
proved removal of cluster galaxies from his background galaxy 
catalogue, increasing both r2oo and c for his best-fit NFW 
model and better agreement of the weak lensing mass profile 
with that derived from strong lensing. We will include colour 
selection of the different galaxy catalogues in a forthcoming 
paper aimed at comparing in great detail all the mass estimates 
at different scales in Abell 1689: velocity distribution of the 
galaxies, X-ray mass maps, strong lensing in the center of the 
cluster and weak lensing at larger scales. Provided the dynam- 
ics of the cluster is well understood this should give a consis- 
tent picture of its mass distribution and components. This is the 
main goal of the pan-chromatic survey which is conducted by 
our group on intermediate redshift X-ray clusters. 

Finally, we will present a global study of our results based 
on the application of the present methodology to the whole 
cluster catalogue, with a discussion of the statistical properties 
of such clusters. A better understanding of the global properties 
of the mass distribution in rich clusters of galaxies will provide 
profound insights into the growth of structure in the Universe. 
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